Parameters of the Effective Singlet-Triplet Model for Band Structure of High-T^ 

Cuprates by Different Approaches 
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The present paper covers the problem of parameters determination for High-T c superconductive 
copper oxides. Different approaches, ab initio LDA and LDA+U calculations and Generalized Tight- 
Binding (GTB) method for strongly correlated electron systems, are used to calculate hopping and 
exchange parameters of the effective singlet-triplet model for Cit02-layer. The resulting parameters 
are in remarkably good agreement with each other and with parameters extracted from experiment. 
This set of parameters is proposed for proper quantitative description of physics of hole doped 
High-T c cuprates in the framework of effective models. 
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I. INTRODUCTION 



High-r c superconductive cuprates (HTSC) belong to 
the class of substances where the strong electron correla- 
tions are important. This circumstance and also the fact 
that these substances have non-trivial phase diagrams 
(see e.g. review 0) lead to the difficulties in description 
of HTSC in framework of first principles (ab initio) meth- 
ods, especially in the low doping region. So at present 
moment the most adequate method of HTSC theoretical 
investigations is the model approach. Effective models of 
HTSC (e.g. t— J model) usually contain free parameters 
that could be fitted to experimental data (comparison of 
calculated and experimental Fermi surfaces, dispersion 
curves, etc.) but as soon as one uses model approach 
the question concerning correctness of these parameters 
arises. One of the possible ways to answer this question 
is to obtain relations between parameters of some effec- 
tive model and microscopic parameters of the underlying 
crystal structure. The underlying crystal structure of 
HTSC could be described either by the 3-band Emery 
model 0, or by the multiband p — d model Q . One 
can compare parameters bound to these models with pa- 
rameters obtained by very different approach, e.g. with 
ab initio calculated parameters. It does not mean that 
the ab initio band structure is correct. Due to the strong 
electron correlations it is for sure incorrect in the low 
doping region where these correlations are most signifi- 
cant. Nevertheless, the single electron parameters are of 
interest and may be compared with the appropriate pa- 
rameters obtained by fitting to the experimental ARPES 
data. 

In the present paper we obtain relations between mi- 
croscopic parameters of the multiband p — d model and 
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parameters of the effective singlet-triplet t — J model for 
hole doped HTSC. Then we compare these parameters 
and the t — J model parameters obtained in the ab initio 
calculations. In Section ^ the details of ab initio cal- 
culation within density functional theory are presented. 
In Section 1TTT1 formulation of the effective singlet-triplet 
model as the low-energy Hamiltonian for the multiband 
p—d model with generalized tight-binding (GTB) method 
applied are described. In both methods the parent in- 
sulating compound La-iCuO^ is investigated. Parame- 
ters are obtained at zero doping because within the GTB 
method the evolution of the band structure with doping 
is described only by changes in the occupation numbers 
of zero-hole, single-hole, and two-hole local terms, while 
all parameters are fitted in undoped case and therefore 
fixed for all doping levels. The resulting parameters of 
both approaches (GTB and ab initio) are in very good 
qualitative and quantitative agreement with each other 
and parameters extracted from experiment. Also, these 
parameters are in reasonable agreement with t — J model 
parameters used in the literature. We conclude that the 
obtained set of model parameters should be used in effec- 
tive models for proper quantitative description of HTSC 
in the whole doping region. 



II. AB INITIO PARAMETERS CALCULATION 

The band structure of LaiCuO^ was obtained in 
frames of the linear muffin-tin orbital method in tight- 
binding approach [f| (TB-LMTO) withing the local den- 
sity approximation (LDA) . The crystal structure data Q 
corresponds to tetragonal La^CuO^. The effective hop- 
ping parameters t p were calculated by the least square fit 
procedure to the bands obtained in LDA calculation 
The effective exchange interaction parameters J p were 
calculated using the formula derived in :9j, where the 
Green function method was used to calculate J p as sec- 
ond derivative of the ground state energy with respect to 
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the magnetic moment rotation angle via eigenvalues and 
eigenf unctions obtained in LDA+C/ calculation The 
LDA+C/ approach allows to obtain the experimental an- 
tiferromagnetic insulating ground state for the undoped 
cuprate in contrast to the LDA approach which gives a 
nonmagnetic metallic ground state |JJj| . The Coulomb 
parameters U = 10 eV and J — 1 eV used in LDA+/7 
calculation were obtained in constrained LSDA supercell 
calculations fll| . 



III. GTB METHOD AND FORMULATION OF 
THE EFFECTIVE SINGLET-TRIPLET MODEL 

The t - J [3 and Hubbard [3 models are widely 
used to investigate HTSC compounds. While using these 
models one, in principle, can catch up qualitatively es- 
sential physics. The parameters in these models (i.e. the 
hopping integral t, antifferomagnetic exchange J, Hub- 
bard repulsion U) are typically extracted from experi- 
mental data. Thus, these parameters do not have a di- 
rect microscopical meaning. More consequent approach 
is to write down the multiband Hamiltonian for the real 
crystal structure (which now includes parameters of this 
real structure) and map this Hamiltonian onto some low- 
energy model (like t — J model) . In this case parameters 
of real structure could be taken from the ab initio calcu- 
lations or fitted to experimental data. 

As the starting model that properly describes crystal 
structure of the cuprates it is convenient to use 3-band 
Emery p—d model 2 , 3] or the multiband p—d model . 
The set of microscopic parameters for the first one was 
calculated in [lj, 03 ■ While this model is simplier then 
the multiband p — d model it lacks for some significant 
features, namely importance of d z 2 orbitals on copper 
and p z orbitals on apical oxygen. Non-zero occupancy of 
d z 2 orbitals pointed out in XAS and EELS experiments 
which shows 2-10% occupancy of d z 2 orbitals |l& ll7| and 
15% doping dependent occupancy of p z orbitals [18[ in 
all HTSC of p-type (hole doped). In order to take into 
account these facts the multi-band p — d model should be 
used: 

Hpd = ^ (e A - /i)«/Aa + Y T /^ C /Aa C <A'<r 

f,X,<7 <f,g>X,X',a 

+ \ Y Y V f X 9 X ' C %<rt C f^s4x'a 2 C 9 yoO,) 

f,g,X,X' ^1,(72,^3,(74 

where Cf\ a is the annihilation operator in Wannier rep- 
resentation of the hole at site / (copper or oxygen) at 
orbital A with spin a, rif\ a — c^ Xa Cf\ a . Indexs A run 
through d x 2_ y 2 = d x and d 3z 2_ r 2 = d z orbitals on 
copper, p x and p y atomic orbitals on the plane oxygen 
sites and p z orbital on the apical oxygen; ex - single- 
electron energy of the atomic orbital A. includes 
matrix elements of hoppings between copper and oxygen 
(tpd for hopping d x <-> p x ,p y ; t pd /\/i for d z <-> p x ,p y ; 
t' pd for d z <-> p z ) and between oxygen and oxygen (t pp 
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FIG. 1: GTB method dispersion (doping concentration 
x = 0) of the top of the valence band and the bottom of 
the conduction band divided by the insulating gap is shown. 
Horizontal dashed lines marks the in-gap states which spectral 
weight is proportional to x. Points with error bars represents 
experimental ARPES data for Sr 2 Cu0 2 Cl 2 



for hopping p x <-> p y ; t' pp for hopping p x ,p y <-► p z ). 

The Coulomb matrix elements includes intraatomic 
Hubbard repulsions of two holes with opposite spins on 
one copper and oxygen orbital (Ud, Up), between different 
orbitals of copper and oxygen (Vd, V p ), Hund exchange 
on copper and oxygen (Jd, J p ) and the nearest-neighbor 
copper-oxygen Coulomb repulsion V p d- 

GTB method 19] consist of exact diagonalization of 
intracell part of p — d Hamiltonian (JIJ and perturbative 
account for the intercell part. For La2- x Sr x CuOi the 
unit cell is CuOq cluster, and a problem of nonorthogo- 
nality of the molecular orbitals of adjacent cells is solved 
by an explicit fashion namely by constructing the rele- 
vant Wannier functions on a five-orbitals initial basis of 
the atomic states [2(], EI • In a new symmetric basis the 
intracell part of the total Hamiltonian is diagonalized, 
allowing to classify all possible effective quasiparticle ex- 
citations in CuO^Tuane according to a symmetry. 

Calculations j2(J, Um °f * ne quasiparticle dispersion 
and spectral intensities in the framework of multiband 
p — d model with use of GTB method are in very good 
agreement with ARPES data on insulating compound 
Sr 2 Cu0 2 Cl 2 HIH (see Fig.©. 

Other significant results of this method are 0> US : 

i) Pinning of Fermi level in Lai- x Sr x CuO^ at low con- 
centrations was obtained in agreement with experiments 
[2(||22]]. This pinning appears due to the in-gap state, 
spectral weight of this state is proportional to doping con- 
centration x and when Fermi level comes to this in-gap 
band then Fermi level "stacks" there. In Fig. [21 the dop- 
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FIG. 2: Dependence of chemical potential shift A/x 
on concentration of doping x for Nd2- x Sr x CuOi and 
La2-xSr x CuOi. Straight lines are results of GTB calcula- 
tions, filled circles with error bars are experimental points 



ing dependence of chemical potential shift Afi for n-type 
High-T c Nd 2 - x Sr x CuO A (NCCO) and p-type High-T c 
La2- x Sr x CuOi (LSCO) is shown. The localized in-gap 
state exist in NCCO also for the same reason as in LSCO, 
but its energy is determined by the extremum of the band 
at (n/2, 7r/2) point and it appears to be above the bottom 
of the conductivity band. Thus, the first doped electron 
goes into the band state at the (it, 0) and the chemical 
potential for the very small concentration merges into the 
band. At higher x it meets the in-gap state with a pin- 
ning at 0.08 < x < 0.18 and then /i again moves into the 
band. The dependence fi(x) for NCCO is quite asym- 
metrical to the LSCO and also agrees with experimental 
data [261 ]. 

ii) Experimentally observed evolution of Fermi 
surface with doping from hole- type (centered at (7T,7r)) 
in the underdoped region to electron-type (centered at 
(0,0)) in the overdoped region is qualitatively reproduced 
in this method. 

iii) Pseudogap feature for Lai- x Sr x CuO ^ is obtained 
as a lowering of density of states between the in-gap state 
and the states at the top of the valence band. 

Above results was obtained with the following set of 
the microscopic parameters: 



u d = v d 



0, EA 



t pp = 0.46, 



9, J d = 1,J P = 0, U p 



= 1.5, 
= 0.58, 
= K = 



the difference ex — es between energy of two-particle sin- 
glet and two-particle triplet depends strongly on various 
model parameters, particularly on distance of the api- 
cal oxygen from the planar oxygen, energy of the api- 
cal oxygen, difference between energy of d 2 2-orbitals and 
dj.2_j,2-orbitals. For the realistic values of model param- 



e p , = 0.45, 
t„ = 0.42, (2) 
V pd 



1.5. 



As the next step we will formulate the effective model. 
Simplest way to do it is to completely neglect contribu- 
tion of two-particle triplet state 3 Bi g . Then there will 
be only one low-energy two-particle state - Zhang-Rice- 
type singlet 1 A\ g - and the effective model will be the 
usual t — J model. But in the multiband p — d model 



eters Et ~ £s is close to 0.5 eV J21J, 32] contrary to the 
3-band model with this value being about 2 eV (this 
case was considered in [2^, H3)- To take into account 
triplet states we will derive the effective Hamiltonian for 
multiband p — d model by exclusion of the intersubband 
hopping between low (LHB) and upper (UHB) Hubbard 
subbands, similar to |12| . 

As the Hubbard X-operators X^ q = \p) (q\ on site / 
represents natural language to describe strongly corre- 
lated electron systems in the rest of the paper we will 
use these operators. The X-operators are constructed in 
the Hilbert space that consists of a vacuum nh = state 
|0), single-hole |er) = {| f),| J.)} state of b\ g symmetry, 
two- hole singlet state \S) of 1 A\ g symmetry and two- 
hole triplet state \TM) (where M = +1,0,-1) of 3 B lg 
symmetry. 

We write the Hamiltonian in the form H = Hq + 
Hi, where the excitations via the charge transfer gap 
Erf are included in H±. Then we define an operator 
H(e) = Hq + eHi and make the unitary transformation 

H(e) = exp(—ieSjH(e)axp(ieS^. Vanishing linear in 
e component of H(e) gives the equation for matrix S: 



Hi 



+ i 



Hq, S 



= 0. The effective Hamiltonian is ob- 



tained in second order in e and at e = 1 is given by: 



H = H + -i 



Hi,S 



(3) 



Thus, for the multiband p—d model in case of electron 
doping (n-type systems) we obtain the usual t— J model: 



5>*. 

f,<> 

<f,g> 



era | \ ^ ,00 y-o-0 Y 0a 
f + 2^ t f9 Ji f A 9 
</,9>,<r 



s f s g 



1 



7 n f n a 



(4) 



where Sf are spin operators and n/ are number of parti- 
cles operators. The Jf g — 2 f*/g) /Ect is the exchange 

integral, E c t is the energy of charge-transfer gap (similar 
to U in the Hubbard model, E ct w 2 eV for cuprates). 
Chemical potential /i is included in e\. 

For p-type systems effective Hamiltonian has the form 
of a singlet-triplet t — J model [3]| : 



H = H + H t + J f9 (SfS g ~ \np 

<f,g> V 



(5) 



where Hq (unperturbated part of the Hamiltonian) and 
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TABLE I: Parameters of the effective singlet-triplet model for 
p-type cuprates obtained in the framework of GTB method 
(all values in eV). 



p 


,00 

L p 


j.3S 
L P 


• OS 

L p 


«TT 
l P 


,ST 
L P 


J P 


(0,1) 


0.373 


0.587 


-0.479 


0.034 


0.156 


0.115 


(1,1) 


0.002 


-0.050 


0.015 


-0.011 





0.0001 


(0,2) 


0.050 


0.090 


-0.068 


0.015 


0.033 


0.0023 


(2,1) 


0.007 


0.001 


-0.006 


-0.004 


0.001 






TABLE II: Comparison of ab initio parameters and param- 
eters, obtained in the framework of GTB method (all values 
in eV). 







b initio 


GTB method 


p 




J P 


t P 


J P 


(0,1) 


0.486 


0.109 


0.587 


0.115 


(1,1) 


-0.086 


0.016 


-0.050 


0.0001 


(0,2) 


-0.006 





0.090 


0.0023 


(2,1) 








0.001 






H t (kinetic part of H) are given by the expressions: 



Ho 



E 



£1 ^ + Z2SX-= S + S2T X 



TMTM 
f 



M 



£ { 

>,<T 

(aV2X 



t fg x f x g 



<f,9>,<? 



,TT 

t fg 



tf g 2a lh 



X 



-x c „ 



Upper indexes of hopping integrals (0,S,T) corresponds 
to excitations which accompanied by hopping from site 
/ to g, i.e. in Hamiltonian one have the following terms: 
The relation between these effec- 



MN vM. 



Xf°X: 



rN 



^ "fg 

<f,g>,v 

tive hoppings and microscopic parameters of multiband 
p — d model is as follows: 



tf g = -2t pd [if g 2uv - 2t pp v lg v 2 



+ SS 

■ ~ 

MS _ 

t fg ~ 

,TT _ 

l fg ~ 



j-ST 



-2t pd ^fg(vj x + U7b) - 2t pp v f gVjb, 

It 

-j=-\f g 2 lalz + 2t pp v fgl 2 a - 2t' pp \ fg 2 lplai 

It 

—^£,fgl z + 2t pp Xfgla ~ 2t' pp £ f g"f p , 



(6) 



The factors n, v, A, £, x are the coefficients of Wannier 
transformation made in GTB method and u, v, j a , jb, 
7zi lp are the matrix elements of annihilation-creation 
operators in the Hubbard X-operators representation. 

The resulting Hamiltonian (JSJ) is the generalization of 
the t — J model to account for two-particle triplet state. 
Significant feature of effective singlet-triplet model is the 
asymmetry for n- and p-type systems which is known 
experimentally. So, we can conclude that for n-type sys- 
tems the usual t — J model takes place while for p-type 
superconductors with complicated structure on the top 
of the valence band the singlet-triplet transitions plays 
an important role. 

Using set of microscopic parameters (j2J) in the Table H] 
we present numerical values of the hopping and exchange 
parameters calculated according to 10 • 



IV. COMPARISON OF PARAMETERS 

The resulting parameters from ab initio and GTB 
calculations are presented in Table [H] Here p is the con- 
necting vector between two copper centers, t p is the hop- 
ping parameter (in the effective singlet-triplet model it 
is equal to t p s , see ©), J p is the antiferromagnetic ex- 
change integral. 

As one can sec, despite slight differences, parameters of 
both methods are very close and show similar dependence 
on distance. It is worth to mention that both methods 
give disproportionality between t p and J p . In the usual 
t — J model proportionality (J p — 2t 2 p /U) takes place 
as soon as this t — J model obtained from the Hubbard 
model with Hubbard repulsion U. In case of the singlet- 
triplet model the intersubband hoppings t® s which deter- 
mines value of J p is different 10 from the intrasubband 
hoppings t p s which determines t p . This fact resulted in 
more complicated relation between t p and J p . 

In the framework of LDA band structure of 
Y Ba2CuOjj rX and withing orbital projection approach 
it was shown [3j| that the 1-band Hamiltonian reduced 
from eight-band Hamiltonian should include not only 
the nearest- neighbor hopping terms (i), but also sec- 
ond (i') and third (t") nearest-neighbors hoppings. In 
GTB method the dependence of hoppings t p on dis- 
tance automatically results from the distance dependence 
of coefficients of Wannier transformation made in this 
method (see Eq. ©). In order to show correspondence 
between results of various authors, we present compar- 
ison of our parameters and parameters, widely used by 
different groups in Table ITTT1 

The parameters extracted from experimental data are 
listed in columns I- VI of Table II I II LDA calculated pa- 
rameters are presented in columns VII and VIII. Our re- 
sults for hoppings are in the best agreement with columns 
III, VII and VIII. This similarity is not surprising. In 
LDA calculations the bandwidth of strongly correlated 
electron systems are usually overestimated due to the 
lack of proper account of strong Coulomb repulsion of 
electrons. But it is well known that the Fermi surface 
obtained by this method is in very good agreement with 
experiments. The main contribution to the shape of the 
Fermi surface comes from kinetic energy of the electrons 
(hopping parameters), so values of hoppings should be 
properly estimated by LDA calculations (columns VII, 
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TABLE III: Comparison of calculated parameters and parameters, used in the literature. 







0" 


P 


IF 


IIP 


IV C 


V c 


VI C 


VlP 


VHP 


IX e 


X e 




LSCO 


LSCO 


LSCO 


LSCO 


LSCO 


Bi2212 


YBCO 


SCOC 


YBCO 


LSCO 


LSCO 


YBCO 


quantity 


here 


here 


[34] 


[35, 36] 


[37, 38] 


[37, 38, 39] 


[40] 


[4JJ 


[33] 


[42] 


[43] 


f43] 


t, eV 


0.587 


0.486 


0.416 


0.35 


0.35 


0.35 


0.40 


0.40 


0.349 


0.43 






t'/t 


-0.085 


-0.18 


-0.350 


-0.20 


-0.12 


-0.34 


-0.42 


-0.35 


-0.028 


-0.17 






t"/t 


0.154 


0.012 




0.15 


0.08 


0.23 


-0.25 


0.25 


0.178 








J, eV 


0.115 


0.109 


0.125 


0.14 


0.14 


0.14 


0.17 


0.12 






0.126 


0.125, 0.150 


J/\t\ 


0.196 


0.224 


0.300 


0.40 


0.40 


0.40 


0.43 


0.30 











a GTB method parameters 

b ab initio parameters obtained in present paper 
c parameters obtained by fitting to experimental data 
d ab initio parameters 

^parameters obtained from two-magnon Raman scattering 



VIII). In works by Tohayama and Maekawa [3?], S3 (col- 
umn III) the parameters was obtained by fitting the 
LSCO tight-binding Fermi surface to the experimental 
one. This procedure should give the same values as an 
LDA calculations and, as one can see, it does. By the 
same technique the parameters for Bi2Sr2CaC'u20s+ x 
(Bi2212, column IV) was obtained [S3,E3| • These param- 
eters arc different from LSCO case and present paper and 
the most straightforward explanation is the more com- 
plicated structure of the Fermi surface of Bi2212 com- 
pound. In the present paper the single-layer (LSCO-like) 
compounds are considered and effects of multiple C11O2- 
planes (i.e. bi-layer splitting) are neglected. The dif- 
ference between our hoppings and hoppings in column 
V appeared due to the same reason (in Ref. 0] the 
Y BCL2CU3O6 insulating compound was investigated). 

In the last two columns of Table ITTT1 the antiferromag- 
netic exchange parameter J obtained from two-magnon 
Raman scattering analysis by moment expansion (LSCO, 
column IX) and spin- wave theory (YBCO, column X) 
are presented (for details see review ji^ and references 
therein). Our values of J (column 0) are in good agree- 
ment with values extracted from experiments and similar 
to listed in columns I- VI. 

In the paper 01 the Hcisenberg Hamiltonian on the 
square lattice with plaquette ring exchange was investi- 
gated. The fitted exchange interactions J = 0.151 eV, 
J' = J" = 0.025J gave values for the spin stiffness 
and the Neel temperature in excellent agreement with 
experimental data for insulating compound La2CuOi. 
In GTB calculations J = 0.115 eV, J' = 0.0009J and 
J" = 0.034 J. The values of J are close to each other but 
different. This difference could be explained by the fact 
that authors of 44] used the Heisenberg Hamiltonian and 
inclusion of hopping term should renormalize presented 
values of exchange interactions. Agreement between J" 
in GTB calculations and Ref. 0] is good but values of 
J' is completely different. The last issue could be ad- 
dressed to oversimplification of calculations in 0] - the 
authors put J' — J" by hand to restrict number of fitting 
parameters. 

Now we will discuss the difference between our param- 




-2 



-2.5 L 1 1 ^ J 

(0,0) (71,71) (7C,0) (0,0)[(7C,0) (0,7C) 

k 

FIG. 3: Dispersion curves on top of the valence band for 
effective singlet-triplet model (singlet subband - solid line, 
triplet subbands - dotted lines) and t — t' — J model (dash- 
dotted line) at optimal doping x = 0.15, dashed line denotes 
self-consistently obtained chemical potential /1. 



eters and parameters in columns I,II,VI and column IV 
(SCOC). The hoppings in mentioned papers were ob- 
tained by fitting t — t' — t" — J model dispersion to 
the experimental ARPES spectra [22, Hij for insulating 
Sr2CuC>2Cl2- We claim that the discrepancy between 
GTB method results and t — t' — t" — J model results 
stems from absence of singlet-triplet hybridization in lat- 
ter model. This statement can be proved by comparison 
of "bare" t — t' — J model Q dispersion and disper- 
sion of singlet-triplet t — t' — J model JSJ. The para- 
magnetic non-superconductive phase was investigated in 
Hubbard-I approximation both in the singlet-triplet and 
t — t' — J models. Results for optimal doping (concentra- 
tion of holes x — 0.15) are presented in Fig. [3] 

There is strong mixture of singlet and triplet bands 
along (0,0) — (7r,7r) and (tt,0) — (0,0) directions due to 
the t ST matrix element (see ©) in both paramegnetic 
(Fig. I3J) and in antiferromagnetic phases (Fig. QJ. It is 
exactly the admixture of the triplet states that deter- 
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mines in our approach the dispersion and the ARPES 
data in the undoped SCOC at the energies 0.3 -j- 0.4 eV 
below the top of the valence band, where the t — t' — J 
model |34| failed, and the additional t" parameter have 
been involved in the t — t' — t" — J model |35|, |3JJ • In our 
approach this parameter is not as necessary as in " bare" 
t — t' — J model since the singlet-triplet hybridization is 
included explicitly. 

In Ref. |45[ the t — t' — t" — J model was also used to 
describe dispersion of insulating SryCuOyClp. and the set 
of parameters was the same as in Refs. [371 133]. But the 
authors of Ref. |4jjj used completely different definition of 
hopping parameters: in their paper the t' term stands for 
hopping between two nearest neighboring oxygens and 
the t" term stands for the hopping between two oxygens 
on the two sides of Cu. Such definition is completely 
different from used in other cited papers where t, t', t" 
terms stands for hoppings between plaquettes centered 
on copper sides and we can't make comparison with their 
results. 

The analysis of Table II I II data gives the following 
ranges for different parameters: 0.350 -j- 0.587 eV for t, 
-0.420 ~ -0.028 for t'/t, 0.012 4- 0.250 for t" /t with the 
exception of value in Ref. [40| and 0.115-^0.150 eV for J. 
In general one can see close similarity in the first neigh- 
bor hopping t and interaction J for the different methods 
and materials, and more discrepancy in such subtle pa- 
rameters as t' and t". 



V. CONCLUSION 

One of the significant results of this paper is the rela- 
tionship @ between microscopic parameters and param- 
eters of the effective singlet-triplet model. Thus the ef- 



fective model parameters are not free any more and have 
a direct physical meaning coming from the dependence 
on microscopic parameters. The parameters of the ef- 
fective singlet-triplet model were obtained both from ah 
initio and model calculations. Model calculations were 
performed in the framework of GTB method for insu- 
lating single-layer copper oxide superconductor. The ah 
initio calculations for LaiCuO^ were done by conven- 
tional LDA TB-LMTO method. The agreement between 
parameters is remarkably good. Obtained parameters 
are also in good agreement with widely used parameters 
of the t — t' — t" — J model but some difference exists. 
This difference is attributed to the neglect of triplet ex- 
citations in simple t — t' — t" — J model. After careful 
analysis we proposed the set of parameters for effective 
models (e.g. t—t'—t"—J model or effective singlet-triplet 
model) for proper quantitative description of physics of 
hole doped High-T c cuprates. 

Acknowledgments 



M.M.K., V.A.G., and S.G.O. thank the Free Univer- 
sity of Berlin for hospitality during their stay. This 
work has been supported by INTAS grant 01-0654, 
Joint Integration Program of Siberian and Ural Branches 
of Russian Academy of Science, Russian Foundation 
for Basic Research grant 03-02-16124, Russian Federal 
Program "Integratsia" grant B0017, Program of the 
Russian Academy of Science "Quantum Macrophysics" , 
and Siberian Branch of Russian Academy of Science 
(Lavrent'yev Contest for Youth Projects), RFFI-01-02- 
17063 (VIA, IAN, MAK), RFFI (MAS)-03-02-06126 
(IAN), Ural Branch of Russian Academy of Science for 
Young Scientists (IAN, ZVP), Grant of the President of 
Russia MK-95.2003.02 (IAN). 



[1] Z.-X. Shen and D.S. Dessau, Phys. Rep. 253, 1 (1995); E. 

Dagotto, Rev. Mod. Phys. 66, 763 (1994); A.P. Kampf, 

Phys. Rep. 249, 219 (1994) 
[2] V.J. Emery, Phys. Rev. Lett. 58, 2794 (1987) 
[3] CM. Varma et al., Solid State Commun. 62, 681 (1987) 
[4] Yu. Gaididei and V. Loktev, Phys. Status Solidi B 147, 

307 (1988) 

[5] O.K. Andersen and O. Jepsen, Phys. Rev. Lett. 53, 2571 
(1984) 

[6] O.K. Andersen, Z. Pawlowska and O. Jepsen, Phys. Rev. 

B 34, 5253 (1986) 
[7] J.D. Axe and M.K. Crawford, J. Low Temp. Phys. 95, 

271 (1994) 

[8] V. I. Anisimov et al., Phys. Rev. 66, 100502 (2002) 
[9] A.I. Lichtenstein et al., J. Magn. Mag. Matter 67, 65 
(1987); A.I. Lichtenstein, V.I. Anisimov and J. Zaanen, 
Phys. Rev. B 52, R5467 (1995) 
[10] V.I. Anisimov, J. Zaanen and O. Andersen, Phys. Rev. B 
44, 943 (1991); V.I. Anisimov et al., J. Phys.: Condens. 
Matter 9, 767 (1997) 
[11] O. Gunnarsson et al., Phys. Rev. B 39, 1708 (1989); V.I. 



Anisimov and O. Gunnarsson, ibid. 43, 7570 (1991) 
[12] KA. Chao, J. Spalek and A.M. Oles, J. Phys. C: Sol. 

Stat. Phys. 10, 271 (1977) 
[13] J.C. Hubbard, Proc. Roy. Soc. A 276, 238 (1963) 
[14] M.S. Hybertsen, M. Schluter and N.E. Christensen, Phys. 

Rev. B 39, 9028 (1989) 
[15] A.K. McMahan, J.F. Annett and R.M. Martin, Phys. 

Rev. B 42, 6268 (1990) 
[16] A. Bianconi et al., Phys. Rev. B 38, 7196 (1988); 
[17] H. Romberg et al., Phys. Rev. B 41, 2609 (1990) 
[18] C.H. Chen et al, Phys. Rev. Lett. 68, 2543 (1992) 
[19] S.G. Ovchinnicov and I.S. Sandalov, Physica C 161, 607 

(1989) 

[20] VA. Gavrichkov et al., Phys. Rev. B 64, 235124 (2001) 
[21] VA. Gavrichkov et al., Zh. Eksp. Teor. Fiz. 118, 422 

(2000); [JETP 91, 369 (2000)] 
[22] B.O. Wells et al., Phys. Rev. Lett. 74, 964 (1995) 
[23] C. Diirr et al., Phys. Rev. B 63, 014505 (2000) 
[24] AA. Borisov, V.A. Gavrichkov and S.G. Ovchinnikov, 

Mod. Phys. Lett. B 17, 479 (2003) 
[25] AA. Borisov, V.A. Gavrichkov and S.G. Ovchinnikov, 



7 



Zh. Eksp. Teor. Fiz. 124, 862 (2003); [JETP 97, 773 
(2003)] 

[26] N.Harima et al, Phys. Rev B 64, 220507(R) (2001) 
[27] A. Ino et al., Phys. Rev. Lett. 79, 2101 (1997) 
[28] A. Ino et al., Phys. Rev. B 65, 094504 (2002) 
[29] J. Zaanen, A.M. Oles, P. Horsch, Phys. Rev. B 46, 5798 
(1992) 

[30] R. Hayn et al., Phys. Rev. B 47, 5253 (1993) 

[31] M. Korshunov and S. Ovchinnikov, Fiz. Tv. Tela 43, 399 

(2001) [Phys. Sol. State 43, 416 (2001)] 
[32] R. Raimondi et al., Phys. Rev. B 53, 8774 (1996) 
[33] O.K. Andersen et al., J. Phys. Chem. Solids 56, 1573 

(1995) 

[34] A. Nazarenko et al., Phys. Rev. B 51, 8676 (1995) 

[35] V.I. Belinicher, A.I. Chernyshev and V.A. Shubin, Phys. 

Rev. B 53, 335 (1996) 
[36] V.I. Belinicher, A.I. Chernyshev and V.A. Shubin, Phys. 



Rev. B 54, 14914 (1996) 
[37] T. Tohayama and S. Maekawa, Supercond. Sci. Technol. 
13, R17 (2000) 

[38] T. Tohayama and S. Maekawa, Phys. Rev. B 67, 092509 
(2003) 

[39] C. Kim et al., Phys.Rev. Lett. 80, 4245 (1998) 

[40] F.P. Onufrieva, V.P. Kushnir and B.P. Toperverg, Phys. 

Rev. B 50, 12935 (1994) 
[41] R. Eder, Y. Ohta, G.A. Sawatzky, Phys. Rev. B 55, 

R3414 (1997) 

[42] E. Pavarini et al., Phys. Rev. Lett. 87, 047003 (2001) 

[43] W. Brenig, Phys. Rep. 251, 153 (1995) 

[44] A. A. Katanin and A.P. Kampf, Phys. Rev. B 66, 

100403(R) (2003) 
[45] T. Xiang and J.M. Wheatley, Phys. Rev. B 54, R12653 

(1996) 



